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ABSTRACT 

The  use  of  boiling  heat  transfer  in  the  liquid  immersion 
cooling  of  electronic  components  has  always  been  hampered  by 
the  excessive  superheat  necessary  to  initiate  nucleation  in 
the  fluorinated  hydrocarbons  used  as  dielectric  cooling 
fluids.  In  an  attempt  to  overcome  some  of  these  difficulties, 
an  experimental  study  of  the  effects  of  nucleate  pooling 
boiling,  on  the  heat  transfer  from  surface  near  the  boiling 
surface  was  conducted.  An  experimental  chamber  was  constructed 
which  had  a  column  of  four  horizontal  wires  spaced  2.5  cm 
vertically  from  each  other.  The  lowest  wire  was  progressively 
heated  from  the  natural  convection  region  through  the  nucleate 
boiling  region.  A  study  was  made  of  the  effects  of  the  boiling 
wake  from  the  lowest  wire  on  heat  transfer  from  the  upper 
wires.  Under  certain  conditions  heat  transfer  enhancements  of 
up  to  30%  were  obtained. 


111 


CI 


TABLE  OF  CONTENTS 

I.  INTRODUCTION    1 

A.  BACKGROUND 1 

B.  PREVIOUS  WORK 3 

C.  OBJECTIVES 7 

II.  EXPERIMENTAL  DESIGN   8 

A.  DESCRIPTION  OF  COMPONENTS   8 

1.  Test  Chamber 8 

2 .  Thermocouples   9 

3.  Heaters   9 

4.  Heat  Exchangers   10 

5.  Aluminum  Cover  Plate  10 

B.  INSTRUMENTATION 10 

1.  Power  to  the  Heaters 11 

2.  Temperature  Measurement   11 

III.  EXPERIMENTAL  PROCEDURE    2  3 

A.  NORMAL  OPERATION  2  3 

B.  DATA  REDUCTION 24 

IV.  RESULTS  AND  DISCUSSION    27 

A.  INDIVIDUALLY  POWERED  27 

B.  PUMPING  EFFECT 2  9 

V.  CONCLUSIONS 50 

VI.  RECOMMENDATIONS 51 

APPENDIX  A.  SAMPLE  CALCULATIONS  52 

APPENDIX  B.  UNCERTAINTY  ANALYSIS   57 

APPENDIX  C.  CALIBRATION  OF  CHROMEL  WIRES   64 

APPENDIX  D.  SOFTWARE   66 

LIST  OF  REFERENCES 7  0 

INITIAL  DISTRIBUTION  LIST    73 


IV 


LIST  OF  FIGURES 

Figure  1.      System  Configuration   13 

Figure  2.      Photograph  of  Test  Apparatus.    14 

Figure  3.  Isometric  View  of  the  Test  Chamber.  .  .  15 
Figure  4.      Cross  Sectional  View  From  the 

Top  of  the  Test  Section 16 

Figure  5.      Schematic    Drawing    of    Test 

Section  in  Vessel 17 

Figure  6.      Voltage     Determination     for 

Thermoelectric  Coolers.    18 

Figure  7.      Technical     Data     for 

Thermoelectric  Coolers.    19 

Figure  8.      Cover  Plate  with  Thermoelectric 

Coolers    20 

Figure  9.      Thermocouple     Connection 

Schematic    21 

Figure  10.     The     Locations     of    the 

Thermocouples  in  the  Chamber.    22 

Figure  11.     Wire   Number    1    Individually 

Powered,  Semi-log 32 

Figure  12.     Number  1  Heat  Flux  vs.  Surface 

Temperature,  Log-log 33 

Figure  13.     Number  3  Shows  Less  Overshoot 34 

Figure  14.     Number  1  at  Different  Times 35 

Figure  15.     All  Wires  Individually  Powered 36 

Figure  16.     Best   Fit   Experimental   Data 

Compared  with  Churchill   (1975) 

Correlation 37 

Figure  17.     Subcooling 38 

Figure  18.     Pumping  Effect  with  Wire  Number 

4  at  42,800  W/m2 39 

Figure  19.     Pumping  Effect  with  Wire  Number 

4  (Detailed  Figure) 40 

Figure  20.     Wire  Number  4  at  Constant  Heat 

Flux  (38,600  W/m2) 41 

Figure  21.     Detail  of  Figure  20 42 


v 


Figure  22.     Wire  Number  3  at  Constant  Heat 

Flux  Value  of  37,490  W/m2 43 

Figure  23.     Wire  Number  3  at  20,430  W/m2 44 

Figure  24a.    Wire  Number  2  at  20,530  W/m2 45 

Figure  24b.    Detail  of  Figure  24a 46 

Figure  25a.    Wire  Number  2  Kept  at  212,700  W/m2.  ...  47 

Figure  25b.    Detail  of  Figure  25a.    48 

Figure  26.     Wire  Number  2  with  Wire  Number  3 

at  Constant  Heat  Flux 49 

Figure  CI.     Chromel  Wire  Calibration 65 


VI 


SYMBOL 

UNITS 

A 

m2 

P 

1/K 

cP 

J/kg 

D 

m 

g 

m/sec2 

Gr 

ND 

h 

W/m2k 

k 

W/mk 

L 

m 

Q 

W 

q 

W/m2 

Pr 

NOMENCLATURE 

DESCRIPTION 

Area 

Thermal  Expansion  Coefficient 

Specific  Heat 

Diameter 

Acceleration  of  gravity 

Grashof  Number 

Heat  Transfer  Coefficient 

Thermal  Conductivity 

Heater  Length 

Heat  Transfer 

Heat  Flux 

Prandtl    Number 

Wall    Superheat,    Tw   -   Tsat 


GREEK  SYMBOLS 


fl  Ns/m2  Dynamic  Viscosity 

v  m2/s  Kinematic  Viscosity 

G  N/m  Surface  Tension 

p  kg/m3  Density 

0)  various  Uncertainty  in  the  Variables 

SUBSCRIPTS 


1 

Liquid 

V 

Vapor 

sat 

Saturation 

sub 

Subcooling 

Vll 


ACKNOWLEDGMENT 

The  author  would  like  to  express  his  gratitude  to 
Professors  M.D.  Kelleher  and  Y.  Joshi  for  their  advice, 
guidance  and  friendship  in  this  endeavor. 

He  wishes  to  thank  the  members  of  the  Mechanical 
Engineering  Support  Shop  for  their  efforts,  especially  to  Jim 
Schof ield. 

Many  thanks  are  due  to  several  other  people  and 
corporations  without  whose  assistance  this  work  could  not  have 
been  completed. 


VI 11 


I .  INTRODUCTION 

A.   BACKGROUND 

Research  on  pool  boiling  in  electronic  cooling  systems  has 
focused  on  three  primary  areas:  (1)  reducing  the  temperature 
excursion  at  incipient  boiling,  (2)  reducing  wall  superheat 
during  nucleate  boiling  and  (3)  enhancing  critical  heat  flux 
(Mudawar,  1990)  .  Since  the  heat  transfer  potential  of  nucleate 
boiling  is  well  known,  this  regime  was  one  natural  selection 
for  examination  as  a  possible  solution  to  this  cooling 
problem.  Nishikawa  (1980)  notes  that  there  are  two  primary 
methods  used  to  promote  nucleate  boiling.  The  first  is  to 
treat  the  surface  in  a  manner  that  reduces  its  wettability. 
The  second,  more  promising  one,  is  to  manufacture  a  surface 
with  cavities  which  trap  and  hold  vapor  keeping  the  nucleation 
site  active. 

In  recent  years  the  drive  toward  higher  heat  flux  and 
power  density  in  electronic  equipment  has  turned  attention  to 
the  use  of  dielectric  liquids.  The  3M  Corporation  has 
commercially  produced  its  "Flourinert"  series  of  high 
dielectric,  inert  electronic  cooling  liquids.  In  contrast  to 
water  and  other  conventional  coolants  these  liquids  possess 
relatively  low  thermal  conductivities  and  boiling  points  and 
small  heats  of  vaporization.  They  demonstrate  high  wettability 


on  most  surfaces.  Their  surface  tensions  are  the  lowest  of  any 
known  liquids,  resulting  in  near-zero  wetting  angles.  This  is 
significant  because  it  eliminates  many  potential  nucleation 
sites.  High  dielectric  strength  and  high  resistivity  are  some 
of  the  characteristics  that  make  these  liquids  well  suited  for 
electronic  applications.  Especially,  FC-72  is  known  as 
chemically  stable  and  leaves  no  residue  during  boiling. 

However,  when  using  dielectric  fluids,  incipience 
superheat  excursion  is  pronounced  which  leads  to  very  high 
wall  temperatures.  A  firm  understanding  of  boiling  incipience 
is  required  in  order  to  implement  effective  cooling.  Boiling 
can  produce  very  large-scale  heat  transfer  coefficients  but  it 
also  introduces  some  problems:  (i)  Boiling  restricts  the 
physical  design  of  the  system,  (ii)  A  high  degree  of  superheat 
may  be  required  if  the  surface  is  very  smooth  in  order  for 
nucleation  to  occur.  Due  to  the  high  local  heat  flux  the 
occurrence  of  periodic  boiling  may  result  in  local  turbulence, 
even  damage  of  electronic  components.  (iii)  Impurities 
concentration  may  result  at  the  nucleation  site. 

Subcooling  offers  many  advantages  over  saturated  boiling 
for  cooling  large  arrays  of  computer  chips  due  to  the 
decreased  diameter  and  rapid  collapse  of  departing  bubbles. 

In  this  study  heat  transfer  enhancement  due  to  bubble 
pumping  in  FC-72  was  investigated.  It  is  hoped  that  they  will 
serve  to  identify  some  trends  and  provide  a  basis  for  further 
research  on  actual  components  and  integrated  circuits. 


B.   PREVIOUS  WORK 

Hsu  (1962)  proposed  a  method  for  determining  the  range  of 
active  nucleation  sites  as  a  function  of  temperature  or  heat 
flux.  This  model  revealed  that  the  maximum  and  minimum  sizes 
of  active  cavities  are  functions  of  subcooling,  system 
pressure,  physical  properties  and  the  thickness  of  the 
superheated  liquid  layer. 

Bergles  and  Chu  (1979)  compared  the  nucleate  pool  boiling 
heat  transfer  characteristics  of  three  copper  Union  Carbide 
High  Flux  test  sections  to  a  plain  copper  tube  in  distilled 
water  and  R-113.  The  experiments  were  conducted  with  three 
different  treatments  to  the  test  surface  prior  to  collecting 
data:  subcooling;  aging  by  preboiling  in  the  pool;  and  heating 
the  surface  in  air  to  remove  all  liquid  in  the  pores.  This 
research  indicated: 


1.  Heat  transfer  coefficients  for  the  porous  surfaces 
were  four  to  ten  times  higher  than  for  the  plain 
surfaces . 

2.  There  was  a  significant  temperature  overshoot 
prior  to  the  initiation  of  boiling  and  it  was  not 
very  sensitive  to  aging  or  the  power  increment 
changes  in  water.  These  overshoots  were  sensitive 
to  aging,  initial  subcooling  and  power  increment 
changes  in  R-113. 


Bergles  et  al.  (1968),  conducted  nucleate  pool  boiling 
studies  in  water,  R-113  and  FC-78,  using  a  304  stainless  steel 
tube.  As  heat  fluxes  were  increased,  the  heat  transfer 
coefficient   followed  the   convective   curve  until   a  high 


superheat  was  attained,  and  then  a  distinct  increase  in  the 
heat  transfer  coefficient  was  noted  as  nucleation  occurred. 
This  was  followed  by  an  abrupt  drop  in  all  temperatures  by  as 
much  as  27°F.  It  was  also  noted  that  the  inception  of 
nucleation  could  be  triggered  by  vibrating  the  test  surface. 
The  conclusions  drawn  were: 

The  phenomenon  of  temperature  overshoot  hysteresis  in 
ordinary  liquids  is  due  to  two  causes: 


1.  The  existence  of  metastable  bubbles  which  are 
triggered  only  at  sufficiently  high  disturbance 
levels, 

2.  The  deactivation  of  larger  cavities  by  displacing 
the  vapor  by  liquid  during  subcooling. 


Many  of  the  experimental  reports  published  prior  to  1986 
are  reviewed  by  Bar-Cohen  and  Simon  (1988).  They  noted 
significant  variations  in  the  incipience  superheat  from  one 
report  to  the  next.  The  authors  discussed  possible  mechanisms 
for  delayed  nucleation.  Recently  non-repeatable  and  unsteady 
behavior  of  boiling  incipience  for  a  highly-wetting  dielectric 
liquid  (R-113)  was  documented  by  You  et    al.     (1989). 

The  bubble  behavior  from  the  nucleate  boiling  regime  to 
the  film  boiling  regime  has  not  been  clarified  yet,  except  for 
liquid  helium.  The  transient  boiling  experiment  of  liquid 
nitrogen  at  stepwise  heat  generation  in  a  platinum  wire  was 
carried  out  by  Kunito  Okuyama  (1989). 


McAdams  et  al .  (1949),  reported  the  effect  of  dissolved 
air  on  forced-convection  subcooled  nucleate  boiling  of  water 
where  the  dissolved  air  reduced  the  heat  flux  at  boiling 
initiation . 

Torikai  et  al.  (1970),  used  platinum  wires  immersed  in 
water  to  measure  dissolved  gas  effect  on  the  incipience  wall 
superheat  under  reduced  pressure  conditions  (0.05  -  1.0  Bar)  . 

Surface  enhancement  is  a  very  popular  tool  for  improving 
pool  boiling  performance.  Marto  and  Lepere  (1982)  found  three 
commercially  available  enhanced  boiling  surfaces  which  shifted 
the  boiling  curve  of  FC-72  toward  lower  superheats  compared  to 
a  plain  copper  surface.  In  a  study  by  Berenson  (1962), 
nucleate  boiling  heat  transfer  coefficients  varied  by  600%  due 
to  variation  in  surface  finish.  He  emphasized  the  importance 
of  the  surface  roughness.  Recent  investigations  on  surface 
effects  in  pool  boiling  by  Chowdhury  and  Winterton  (1985)  show 
that  heat  transfer  appears  to  be  unaffected  by  contact  angle. 

The  effects  of  flow  velocity  and  subcooling  on  critical 
heat  flux  from  a  simulated  electronic  chip  attached  to  the 
wall  of  a  vertical  rectangular  channel  was  determined  by 
Mudawar  and  Madox  (1988) .  They  concluded  that  heat  flux  was 
significantly  enhanced  by  increasing  the  degree  of  subcooling. 

The  solubility  of  gas  in  highly-wetting  dielectric  fluids 
is  of  special  concern.  These  fluids  can  contain  up  to  25  times 
more  dissolved  air  by  volume  than  can  water:  48%  air  by  volume 


can  be  dissolved  in  FC-72  compared  to  1.9%  air  by  volume  in 
water  (3M  Product  Manual,  1987). 

Experimental  results  presented  by  You,  Simon  and  Cohen 

(1990)  indicate  that  dissolved  air  plays  an  important  role  at 

incipience  only  when  a  very  high  dissolved  air  content  is 

maintained.  The  surface  tension  value  for  FC-72  was  estimated 

to  have  a  2%  reduction  due  to  high  dissolved  air  content. 

S.M.  You  et  al.  (1989)  conducted  an  experiment  in 
saturated  R-113  fluid  under  atmospheric  pressure  with  two 
types  of  wires.  One  is  0.13  mm  diameter  chromel  wire  (90% 
nickel)  and  the  second  is  a  thin  platinum  film  on  a  quartz 
rod.  They  examined  many  pits  of  various  sizes  using  a  scanning 
electron  microscope  (SEM) . 

Bar-Cohen  (1990)  conducted  experiments  providing  pool 
boiling  characteristics  for  a  number  of  commercially  available 
inert  liquids  (FC-43,  FC-72,  FC-75,  FC-77,  FC-84,  FC-87) .  Data 
were  obtained  with  a  classical  0.25  mm  platinum  wire  heating 
element.  He  concluded  that: 


1.  Nucleate  boiling  characteristics  of  these  fluids 
depend  on  the  heating  surface, 

2.  Saturated   pool   boiling   curves   were   almost 
identical  with  one  another  on  a  particular  heater, 

3.  Superheat  excursions  at  boiling  incipience  appear 
to  always  occur,  but  their  extent  can  vary  widely. 


C.   OBJECTIVES 

The  objective  of  this  work  was  to  investigate  the  effects 
of  the  boiling  wake  from  a  heated  wire  on  the  heat  transfer 
from  wires  placed  within  the  wake.  This  work  was  carried  out 
using  the  dielectric  liquid  FC-72. 


II.  EXPERIMENTAL  DESIGN 

A.   DESCRIPTION  OF  COMPONENTS 

Figure  (1)  shows  the  system  configuration  and  identifies 
all  major  components.  Figure  (2)  is  a  photograph  of  the 
apparatus  . 

The  chamber  consisted  of  a  plexiglass  box  with  an  aluminum 
coverplate  fitted  with  a  rubber  O-ring.  This  vessel  has  side 
wall  heaters  to  allow  the  liquid  to  be  preboiled  for  degassing 
before  the  run,  and  to  maintain  the  liquid  at  saturation 
temperature  throughout  the  run.  The  vapor  was  condensed  and 
returned  to  the  chamber  by  gravity  from  the  aluminum  cover 
plate  condenser  equipped  with  thermoelectric  coolers.  The 
following  is  a  description  of  the  test  apparatus. 

1 .  Test  Chamber 

The  test  chamber  was  constructed  of  plexi-glass  of  1/2 
inch  thickness  and  had  inside  dimensions  of  6  inches  height, 
6  inches  length  and  2  inches  width.  An  isometric  view  of  the 
chamber  is  presented  in  Figure  (3) .  Side  walls  are  glued  to 
each  other. 

It  has  a  1  inch  diameter  hole  and  1  inch  protruding 
pipe  and  a  copper  elbow  tube  in  the  back  face  which  is  the 
pathway  for  the  power  and  thermocouple  wires.  The  chamber  also 
allows   the   test   section  to   slide   in   and  out   for  the 


experiment.  Cross-sectional  and  schematic  drawings  of  the  test 
section  in  the  chamber  are  presented  in  Figures  (4)  and  (5) 
respectively.  The  chamber  also  contained  two  small  (1/4  in) 
pressure  equalization  outlets,  located  at  the  top-right  and 
the  bottom-left  of  the  front  vertical  surface.  It  is  a  simply 
a  closed  tank  with  an  aluminum  cover  plate. 

The  test  chamber  is  attached  to  a  wood  base  in  order 
to  give  better  arrangement  for  cables  and  connections. 

2 .  Thermocouples 

Thermocouples,  0.01  in  (0.254  mm)  in  diameter,  were 
used  to  take  measurements.  The  thermocouples  for  the  plate 
condensers  (thermoelectric  coolers  which  are  located  on  top  of 
the  aluminum  plate)  were  bonded  in  grooves  by  placing  a  small 
drop  of  Eastman  910  adhesive  on  the  thermocouple  bead.  The  two 
thermocouples  then  were  placed  in  the  designated  line  groove. 
Using  a  weight,  pressure  was  applied  until  the  adhesive  set. 
After  allowing  for  three  hours  cure  time,  the  remainder  of  the 
groove  was  filled  with  Omega  Bond  101  Epoxy  and  smoothed  to 
the  level  of  the  aluminum  plate  inside  face.  The  Omega  Bond 
101  Epoxy  experienced  little  shrink  after  curing. 

3 .  Heaters 

The  side  wall  heaters  were  glued  on  either  side  walls 
of  the  chamber.  The  actual  size  of  the  heaters  were  125  x  11 
mm  and  they  were  subjected  to  5  V,  2.5  amperes.  The  power  load 
depended  on  the  combination  of  the  wires  being  heated. 


4 .  Heat  Exchangers 

The  three  thermoelectric  heat  exchangers  measure  1.56 
x  1.56  inches  and  were  placed  on  top  of  the  aluminum  cover 
plate.  The  power  supply  for  these  were  determined  by  the 
related  curves.  The  property  curves  and  technical  data  are 
presented  in  Figures  (6)  and  (7)  respectively. 

5.  Aluminum  Cover  Plate 

An  aluminum  plate  of  1/2  inch  thickness  was  used  as  a 
cover  plate  for  the  test  chamber.  The  inside  face  of  the  plate 
had  two  thermocouple  beads  to  monitor  the  condenser 
temperature.  Three  thermoelectric  condensers  were  placed  on 
top  of  the  aluminum  plate.  A  thin  groove  was  cut  around  the 
edges  of  the  inside  face  for  a  O-ring  gasket.  The  plate  was 
screwed  over  the  test  chamber  with  eight  stainless  steel 
screws  of  1/16  inch  diameter  and  compressed  the  O-ring  to  form 
a  watertight  seal.  Figure  (8)  shows  the  cover  from  the  top 
with  the  condensers  on  it. 

B .   INSTRUMENTATION 

Separate  voltage  taps  across  the  heating  element  and 
across  the  precision  resistor  are  used  to  compute  heater 
resistance  and  heating  element  power.  The  heating  element 
temperature  is  computed  from  its  resistance  via  temperature 
calibration  curve  obtained  prior  to  the  experiment. 
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1 .  Power  to  the  Heaters 

Each  heater  was  run  in  series  with  a  precision 
resistor  that  was  measured  to  have  a  0.1  ohm  resistance.  All 
four  heaters  were  in  parallel  with  a  Hewlett  Packard  model 
number  6289A,  0-40  volts,  0-1.5  amperes,  DC  power  supply.  Both 
the  voltage  across  the  precision  resistor  and  the  chromel  wire 
were  measured  independently.  The  current  to  each  heater  was 
calculated  by  the  product  of  the  heater  voltage  and  heater 
current.  Both  voltages  were  measured  by  a  Hewlett  Packard 
model  3852S  data  acquisition  system  containing  a  Hewlett 
Packard  model  44701A  digital  voltmeter,  all  controlled  by  a 
Hewlett  Packard  model  300  computer. 

2 .  Temperature  Measurement 

The  thermocouples  for.  the  bulk  fluid,  ambient,  ice 
reference  point  and  coolers  inside  plate  temperatures  were 
referenced  individually  to  an  electronic  Ice  Point  Reference 
as  seen  in  Figure  (9) .  Each  reference  thermocouple  was 
connected  such  that  its  constantan  lead  was  connected  to  a 
constantan  lead  of  the  measurement  thermocouple.  The  copper 
leads  of  each  measurement  and  reference  thermocouple  were 
connected  directly  to  a  Hewlett  Packard  model  44705A  20 
channel  relay  multiplexer  and  inserted  into  the  data 
acquisition  system.  The  data  acquisition  system  then  measures 
the  ice  referenced  voltage  of  each  thermocouple.  The  voltages 
were  converted  directly  into  temperatures  in  the  controlling 
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computer  program  by  using  a  fourth  order  polynomial,  fit  to 
the  thermocouple  manufacturer' s  calibration  data  for  copper- 
constantan  thermocouples  (Hazard,  1986) .  The  locations  of  the 
thermocouples  in  the  chamber  are  presented  in  Figure  (10) . 
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Figure  2 .  Photograph  of  Test  Apparatus 
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Figure  3.  Isometric  View  of  the  Test  Chamber 
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Key 


1 .  Washers 

2.  Chromel  Wire 

3.  Screw 

4.  Nut 

5.  Power  Line 


Figure  5.  Schematic  Drawing  of  Test  Section  in  Vessel 


17 


IfldNI   S110A 
5       5!        °o       *       o 


"     D-  Qldmd   1V3H   illVM 
5  * 


LO 


^S         2i  <ft         5-         Q 

iHdNI      S110A 


T 


I       §        3       5 
&"^-  Q3d^nd  1V3H  'MH/nig 


18 


J 


U 
10 


CD 
U 
<D 
H 
0 
O 

o 

•H 

u 

o 
<D 

H 
<D 
O 


<D 

E-" 

U 

0 


-o 

0) 


(TJ 

Q 


T3 
C 
10 


n 


Bl     B      1) 

Xi  "o 

(0    — i 

c  ci 
o 


*-.        o 


w 

.    CD    CO 
.    w   — J 

£  a  ° 


a-g 


CD 


52  £ 


O  V4-1 

•4-"  U 

1-  <" 

«  a 

a  ol 


■2  o 

5  E 

(0        Q) 

££  cx~|* 


>  ^ 

«  -J 
H 
O 
tn    co 

r-l    e. 

H 

a 

u  C 
3  C 
05    t-l 


o  -a 
X    C 


o 

•H 

c 
u 

d) 


a) 
h 

3 
Cn 

■H 

fa 


19 


Ul 
V 

z 
ttf 

< 

® 

<S 

,i 

no 

VI 

2 
o 

2 
ui 

r 

a 

j 
j 
< 

$ 

L9 

§ 

t\l 

1 

L 

1 
T 

<S 

— =a 

£-J 

r 

— 

CO 

1 

■ 

sO 

55- 

"*             r*                      "1 

\ 

0) 

w 

(D 
H 
0 
0 
O 

o 

•H 
U 

U 
0 

H 
® 

i 

-C 

•H 

5* 

(D 
■P 
(Q 
H 

04 

H 
<D 

g 


CO 

0 

u 

& 
•H 
fa 


20 


or 


5° 

c2 

\ 

3i 

end 

<    MJ 

LU  — 

St- 

LU 

& 

—  UJ 

£°i 

uJ 

Uju- 

(J  Ui 

v-a: 

z 

o . 

—  s 

DATA 
UISIT 
YSTE 

O  en 

u 

< 

c: 


LU 

u 

LU 

£T 
LU 

L_ 
LU 

C 

LU 

u 


o 

•H 
4J 

m 
I 

<D 

o 

CO 

c 

0 
•H 

-p 
u 

<D 

c 
c 

0 

u 
<d 

H 

a 

3 
0 

o 

0 


0 

XI 


■H 


21 


COOLERS 


r~ 
h 


».    .<     .    .' 


VAPOR 


7v       S*\ 


0 ! o 

TEST    SECTION/ 
o a 


Figure  10.  The  Locations  of  the  Thermocouples  in  the  Chamber 
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III.  EXPERIMENTAL  PROCEDURE 

A.   NORMAL  OPERATION 

In  carrying  out  the  experiments,  the  procedure  was  as 
follows : 

1.  The  ice  point  reference  device  was  turned  on.  A  mercury 
thermometer  inserted  into  the  ice  base  ensured  that  the  ice 
bath  is  at  0  ±  0.05°C  within  one  hour. 

2.  The  data  acquisition  system  was  then  turned  on.  The 
internal  voltmeter  in  the  data  acquisition  system  requires  one 
hour  warm-up  period. 

3.  The  computer  and  printer  were  turned  on.  The  data 
acquisition  program  contained -in  Appendix  D  was  loaded. 

4.  The  HP  6289A,  0-40  volt,  DC  voltage  supply  was  set  to 
zero  and  energized.  The  current  limiter  on  the  40  volt  power 
supply  was  set  to  the  maximum  current  level. 

A  sampling  of  all  thermocouple  temperatures  was  performed 
and  printed.  The  power  was  increased  incrementally  and  the 
output  from  the  power  calculation  program  was  checked. 

The  bulk  fluid  temperature  as  given  by  the  wire  resistance 
must  match  the  value  given  by  the  bulk  fluid  thermocouples. 
When  the  results  showed  that  the  temperature  stratification  to 
be  less  than  0.10°C  the  experiment  was  ready  to  begin.  The 
power  to  the  test  section  is  then  increased  by  controlling  the 
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test  section  heater  voltage.  The  voltage  was  increased  in  0.2 
volt  increments  until  the  initiation  of  boiling  occurred,  at 
which  time  0.1  volt  increments  were  made  until  the  maximum 
attainable  heat  flux  was  reached.  This  occurred  at  5.7  volt 
for  the  power  supply,  making  use  of  the  HP  6289A,  0-40  volt, 
0-1.5  Amperes  DC  voltage  supply  up  to  8.1  volts,  which 
corresponds  to  the  maximum  value.  The  power  supply  has  a  reset 
function  to  protect  the  device  itself.  The  power  is  then 
decreased  in  the  same  fashion.  At  each  power  setting  the 
system  was  permitted  to  stabilize  for  a  few  minutes  and  the 
following  data  were  recorded:  The  heater  voltage,  precision 
resistor  voltage,  two  thermocouple  indications  of 
thermoelectric  cooler  temperatures,  and  three  thermocouple 
readings  of  pool  temperatures.  This  procedure  was  repeated  for 
every  set  up. 

B.   DATA  REDUCTION 

The  test  wire  temperature  was  found  using  the  calibration 
formula  between  thermal  resistance  and  temperature  value  of 
the  chromel  wire. 

R  =  R0(l+eAT) 
where, 

AT    =       Ts   -    Tb 

Ts   is  the  wall  temperature 

Tb   is  the  bulk  liquid  temperature 

£  is  the  slope  of  the  calibration  curve. 
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Details  of  the  calibration  slope  value  are  presented  in 
Appendix  C.  This  study  utilizes  the  value  of  0.0004  C/ohm 
offered  by  Kaye  (1976) . 

The  starting  value  of  the  wire  resistance  is  accepted  as 
the  R0   in  the  previous  formula. 

The  heat  flux  is  calculated  from  the  following  formula: 

Flux  =  Power /A 
where,  A   =  n   D   L   is  the  total  wetted  surface  area,  and  D   and 
L   are  diameter  and  length  of  the  wire  respectively. 

The  bulk  temperature  of  the  liquid  was  measured  by  taking 
the  arithmetic  average  of  the  three  thermocouple  readings. 

Tbulk  =  (  Tx   +  T2  +  T3  )  /  3 
The  Rayleigh  number  based  on  diameter  is  given  as: 

g$  (Ts-Tb)  d3 

D  va 

where, 


e  = 


vpfCp) 

All  the  properties  in  the  above  equations  are  calculated 
at  the  bulk  temperature. 

Determination  of  Nusselt  number  is: 

d  k 

where  h  is  the  heat  transfer  coefficient  and  was  determined 

from 
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h- * 

I  Ts  -  Tbulk  ) 
where, 

Ts     =       wire  surface  temperature  (°C) 

h       =       heat  transfer  coefficient  (W/m2  •  C) 

q"     =   heat  flux  (W/m2) 

Churchill   and  Chu   (1975)   have   recommended  a  single 

correlation  formula  for  the  natural  convection  from  long 

horizontal  cylinders. 

-   ,  0.381Rao/6  ., 

Nud  =  {  0  .  60  + _,  .  ...  )2 

d  [  1  +  (  0.559/Pr)9/ie  ]8/2/ 

It  is  valid  for  a  wide  Rayleigh  number  range 

10-5  <  Rad  <  1012 

where  the  Prandtl  number  (Pr)  is  defined  as 

a 
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IV.  RESULTS  AND  DISCUSSION 

Data  were  taken  and  the  related  graphs  can  be  grouped  for 
discussion  purposes  as: 

1.  individually  powered  wires 

2.  pumping  effect  runs 

A.   INDIVIDUALLY  POWERED 

Before  investigating  the  pumping  effect  of  the  bubbles, 
individual  boiling  experiments  on  the  wire  heaters  were 
conducted. 

Figure  (11)  shows  wire  1  (the  lowest  wire)  where  distinct 
single  phase  and  nucleate  boiling  regimes  are  observed.  When 
the  heat  flux  was  increased  in  small  steps,  the  single  phase 
curve  reaches  17°C  superheat.  At  about  2  0°C  an  abrupt 
temperature  decrease  occurs  as  nucleate  boiling  begins. 

Figures  (11)  and  (12)  are  essentially  the  same.  The 
coordinates  are  heat  flux  versus  surface  temperature  of  the 
chromel  wire.  Figure  (12)  is  on  the  log-log  scale  while  Figure 
(11)  is  semilog  graph. 

Figure  (13)  is  the  single  run  boiling  curve  of  heater 
number  3.  It  was  the  only  one  to  show  overshoot  of  less  that 
7°C. 
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Figure  (14)  for  wire  number  1  shows  the  same  type  of 
behavior  as  Figure  (11)  .  In  this  case  data  were  taken  at 
different  times  in  order  to  demonstrate  the  repeatability  of 
the  experiments.  Data  points  for  two  different  runs  did  not 
coincide,  but  that  is  because  the  bulk  fluid  temperature  is 
not  the  same  for  all  experiments,  but  depends  on  the  ambient 
temperature  of  the  room. 

Figure  (15)  indicates  all  individually  powered  heaters  on 
the  same  graph.  The  difference  in  the  overshoot  may  be  because 
of  the  surface  geometry  of  heater  number  3. 

The  large  aspect  ratio  of  the  heaters  resulted  in  minimal 
conduction  losses  to  the  supports  and  yielded  a  uniform 
temperature  along  the  length. 

The  boiling  curves  for  all  four  wires  are  practically 
superimposable,  even  at  the  departure  from  nucleate  boiling. 
They  differ  only  in  the  heat  flux  at  which  boiling  incipience 
occurred  and  at  which  full  nucleate  boiling  began. 

Figure  (16)  shows  the  best  fit  to  the  natural  convection 
portion  of  the  experimental  data  for  wire  number  1,  and  is 
compared  with  Churchill  (1975)  correlation  values.  There  is  a 
constant  offset  from  the  correlation  values  because  of  the 
wall  efffect  except  for  the  low  Rayleigh  numbers  of  that  run 
which  also  shows  considerable  scatter  in  the  data. 

The  saturation  temperature  of  the  FC-72  is  about  56°C. 
Tsat  ~  Tbuik  is  the  subcooling.  Figure  (17)  shows  the  effects  of 
subcooling.  At  boiling  incipience,  a  liguid  mass  near  the 
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heated  surface  is  slightly  superheated,  even  though  the  bulk 
fluid  is  subcooled.  Each  run  began  with  zero  heating.  The  heat 
flux  value  was  increased  incrementally,  continuing  up  the 
single-phase  natural  convection  line  until  boiling  was 
initiated.  The  figure  indicates  that  as  the  subcooling 
increases  the  curve  shifts  to  the  left,  displaying  significant 
increases  in  the  overshoot  and  heat  fluxes.  With  bulk  liquid 
at  15°C,  it  shows  large  amount  of  superheat,  probably  a  bubble 
begins  to  form  on  the  wire  and  as  soon  as  the  bubble  grows  a 
bit  it  encounters  subcooled  liquid  and  collapses.  The  liquid 
is  so  cold  that  it  doesn't  allow  the  bubble  to  grow.  As  the 
bulk  temperature  of  the  liquid  gets  warmer,  it  allows  the 
bubbles  to  grow.  Subcooling  results  in  smaller  bubbles,  which 
are  less  likely  to  spread  to  neighboring  nucleation  sites  than 
the  larger  bubbles  associated  with  saturated  boiling. 

B.   PUMPING  EFFECT 

For  this  type  of  experiment,  heater  number  1  was  used  as 
the  pumper. 

Wire  number  4  was  fixed  at  various  heat  flux  values,  while 
heater  number  1  was  increased  in  a  stepwise  manner. 

Figures  (18)  and  (19)  show  the  pumping  effect  of  the 
bubbles  generated  by  wire  number  1  on  the  temperature  of  wire 
number  4.  Clearly  Figure  (19)  gives  better  understanding  by 
way  of  marked  corresponding  points.  As  the  temperature  comes 
closer  to  the  temperature  at  which  the  boiling  incipience 
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starts,  wire  number  4  is  little  effected.  Just  after  wire 
number  1  is  subjected  to  overshoot  it  also  has  similar  a 
experience,  i.e.,  a  sharp  jump  .  Actually,  it  was  the  general 
pattern  for  most  of  the  other  runs.  This  suggests  that  boiling 
can  be  started  by  the  pumping  effects  of  the  rising  bubbles 
which  trigger  the  heater  located  above  it.  The  characteristics 
of  the  bubbles  leaving  the  lower  heater  are  uncertain.  They 
may  be  the  wake  of  the  buoyancy-induced  motion,  especially 
when  the  upper  wire  was  fixed  at  a  high  heat  value.  A  similar 
trend  can  be  observed  from  Figure  (20)  .  One  thing  to  notice  is 
that  heater  number  4  is  fixed  at  a  heat  flux  that  corresponds 
to  a  starting  value  of  surface  temperature  of  50°C  for  this 
run.  This  is  below  the  temperature  that  meets  the  overshoot. 
Less  jump  in  the  temperature  is  associated  with  the  lower 
fixed  heat  flux  value.  If  we  compare  it  with  Figure  (19), 
which  was  at  higher  heat  value  than  this  one,  Figure  (21) 
provides  a  detailed  look  at  Figure  (20) . 

Figure  (22)  shows  wire  number  3  fixed  at  37,490  W/m2.  It 
has  the  points  marked  in  chronological  order.  The  3rd  wire 
heater  is  located  closer  to  the  first  wire  than  is  the  fourth 
wire,  so  it  allows  most  of  the  bubbles  to  reach  the  upper 
heater  with  collapsing  due  to  surrounding  subcooled  liquid  FC- 
72. 

Figure  (23)  shows  another  combination  of  heater  number  3 
with  heater  number  1.  Heater  number  3  is  kept  at  20,450  W/m2. 
This  was  very  low  heat  flux  value  that  the  triggering  effect 
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did  not  show  itself.  Bubbles  still  caused  some  decrease  in  the 
temperature  of  the  upper  wire. 

Figures  (24a)  and  (24b)  show  very  similar  results.  In  this 
case  wire  number  2  is  kept  at  20,530  W/m2.  Since  wire  number 
3  is  the  closest  to  wire  number  1.  Being  so  close  helped  it  to 
lose  some  heat  corresponding  to  a  decrease  in  its  temperature. 

Figures  (25a)  and  (25b)  show  no  effect  from  the  bubbles 
from  below.  Here  heater  number  2  is  kept  at  a  very  high  heat 
flux  of  212,700  W/m2,  which  is  already  boiling.  It  is  not 
sensitive  any  more  to  the  bubbles  produced  by  heater  number  1. 
Consequently  no  jumps  or  apparent  decrease  on  the  surface 
temperature  of  the  number  4  wire  are  seen. 

Figure  (26)  shows  similar  trend  as  the  early  graphs.  In 
this  case  the  bubbles  are  produced  by  wire  number  2  and  the 
effects  on  wire  number  3  are  presented.  Heater  number  3  had 
the  temperature  decrease  as  the  pumper  had  overshoot 
experience . 
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WIRE  #4  AT  CONS.HEAT  FLUX 
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Figure  18.  Pumping  Effect  with  Wire  Number  4  at  42,800  W/m2 
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Figure    19.    Pumping  Effect  with  Wire  Number  4 
(Detailed  Figure) . 
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WIRE  #4  AT  CONST.  HEAT  FLUX 
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Figure  20.  Wire  Number  4  at  Constant  Heat  Flux  (38,600  W/m2) 
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Figure  21.  Detail  of  Figure  20 
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Figure  22 .  Wire  Number  3  at  Constant  Heat  Flux  Value  of 

37,490  W/m2. 
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WIRE  #3  AT  CONS.  HEAT  FLUX 
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Figure   23.    Wire  Number  3   at   20,430  W/m2. 
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WIRE  #2  AT  CONS.  HEAT  FLUX 
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Figure   24a.    Wire   Number   2   at   20,530   W/m2 
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WIRE  #2  AT  CONS.  HEAT  FLUX 
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WIRE  #2  AT  CONS.HEAT  FLUX 
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Figure  25a.  Wire  Number  2  Kept  at  212,700  W/m2 
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WIRE  #2  AT  CONS.HEAT  FLUX 
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Figure   25b.    Detail   of  Figure   25a 
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WIRE  #3  AT  CONS.  HEAT  FLUX 
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V.  CONCLUSIONS 

The  figures  give  the  idea  of  the  pumping  effects.  When  the 
pumping  heater  is  powered  progressively,  the  upper  heater  will 
react  as  indicated  in  the  following. 

If  the  upper  wire  is  located  close  to  the  lower 
pumper,  it  will  meet  most  of  the  bubbles.  It  will  lose  some 
heat  depending  on  how  close  it  is. 

•  The  initial  heat  flux  value  of  the  upper  heater  is 
also  important.  The  amount  of  temperature  jump  which 
simultaneously  occurs  when  the  pumper  has  overshoot  depends  on 
how  warm  the  upper  wire  is. 

•  If  the  upper  heater  already  has  bubbles,  this  means  it 
had  incipience  before,  and  the  bubbles  will  not  help  it  lose 
more  energy. 


50 


VI .  RECOMMENDATIONS 

1.  Use  platinum  wires  rather  than  chromel  wires  to  get 
better  temperature  measurement. 

2.  Before  the  start  of  each  run,  boil  the  pool  for  a 
sufficient  time  to  degas  the  liquid. 

3.  Modify  the  software  for  taking  data  by  adding  a  WAIT 
command  to  obtain  the  time  function. 

4.  An  addition  of  a  foil  type  heater  placed  horizontally 
at  the  bottom  of  the  chamber  may  help  create  a  uniform 
temperature  throughout  the  bulk  fluid. 

5.  An  extension  of  tygon  tubing  of  about  one  foot  length 
would  take  care  of  overflow  and  hydrostatic  head  determined  by 
the  liquid  level  inside  the  filling  tank. 

6.  The  present  study  deals  with  only  subcooling  effects. 
Saturated  pool  boiling  should  also  be  studied.  In  order  to 
have  constant  temperature  in  the  bulk  fluid,  the  thermo- 
electric condensers  and  the  bulk  heaters  can  be  run 
simultaneously . 
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APPENDIX  A.  SAMPLE  CALCULATIONS 

A.   DETERMINATION  OF  INPUT  POWER 

For  the  sample  calculations,  the  output  values  of  wire 
number  1  were  used: 

Power ( 1  )  =  voI5  ( 12)  volt  ( 11 ) 
DTI 

where  volt  (12)  and  volt  (11)  are  the  voltages  across  the 
precision  resistor  and  wire  respectively. 


=  (  °-0048  j( 1.988  )  =  0.96  Watts 


B.  CONVERSION   FROM   VOLTAGES   TO   TEMPERATURES    (°C)    FOR 
THERMOCOUPLES 

The  HP  3852  A  Data  Acquisition  System  channels  1-9  record 
the  voltage  readings  in  millivolts.  The  following  polynomial 
formula  is  for  copper  constantan  thermocouples. 

T  (  °C)  =0.1008  60  91+  (25727.9*  emf )  -   (7  67  343.5*  emf2  ) 

+  (7802.55  96*  emf3  )-  (924748  6589. 6*  emf4)  + 

(  6.  9BE11   *  emfb  )  -  (  2.66E13  *  emf6  )  +   (  3.94E14  *  emf1  ) 

Using  thermocouple  1 

emf   =  0.000779  volt         T=19.72°C 

C.  CALCULATION  OF  WIRE  TEMPERATURE 

The  calibration  curve  for  the  chromel  wire  (Kaye,  1976) 
yields  the  slope  £  of  4.10"4  ohm/°C. 
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R  =  R0[  I  +ET-T0] 
Here  R0   and  R   are  the  resistances  of  the  wire  at  the  reference 

point  and  at  any  point  to  be  measured  .R0  and  ohm(l)  are  the 

corresponding  variables  used  in  the  software. 

Ts(l)    =    ((  (ohm(l)  /RJ    -    1)    /  Kaye    )     +    T0 

where, 

Ts(l)       =       wire  temperature 

Kaye        =       0.0004    ohm/°C    constant    slope 

T0  =       wire    reference    temperature    (°C)  . 

Ts(l)    =    ((  (4 .1009/4 .075)     -    1)     /    4E-4)     +    19.64 

=   35.5°C 

D.  CALCULATION  OF  THE  AVERAGE  BULK  TEMPERATURE 

The  average  temperature  of  the  liquid  was  simply  the 
arithmetic  average  of  the  three  thermocouple  values  immersed 
in  the  bulk  fluid. 

Tbulk  -  (T(l)     +    T(2)     +    T(3))    /   3 

=    (19.9    +    19.54    +    19.6)    /   3   =   19.68°C 

E.  DETERMINATION  OF  HEAT  FLUX 

HFLUX (  1  )  = POWer < 1  ) 


AREA 


where, 


AREA         =       K    *    D    *    L    =     (3.1415926)    D    L 

D,    L        =       Diameter  and  length  of  the  wire. 

AREA        =        (3.1415926)     (0.005)     (0.02542)     (3)    =   3.04E-5   m2 

HFLUX =  ^2;355  =774  60  W /  m2 
3.04E-5 
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F.   FLUID  PROPERTIES  DETERMINATION 

1.  Thermal  Conductivity  k(W/m°C) 

From  Figure  5  of  the  3M  Corporation  Flourinert  Product 
Manual,  the  thermal  conductivity  coefficient  curves  have  been 
determined  to  be 

K=   (  0.6033  -  0.00115  *  Th)  I  10     lmW  =0.1  JL 

b  Cm°C  mK 

at  Tb   =  50°C   it  yields  0.0545  W/m. 

2 .  Liquid  Density  p  (kg/m3) 

Using  the  expression  in  Table  4B  and  constants 
presented  in  Table  AC   of  the  Product  Manual  yields 

p  =  (1.740    -    0.00261    Tb)     *    1000 
Tb  film  temperature  must  be  in  units  of  °C.  At  Tb  =  50°C  this 
yields  p=  1609.5  kg/m3. 

3.  Kinematic  Viscosity  v  (m2/s) 

From  Figure  (3)  and  determining  a  4th  order  curve  fit 
yields 

V  =  (0.54    -    0.0116    *    Tb    +    0.0002085    *    Tb2)     *   IE- 6 
at  Tb   =  50°C  this  yields  V  =  4.8E-7  m2/s. 

4 .  Density  of  Vapor  pv  (kg/m3) 

From  Table  13A  of  the  Product  Manual,  values  for  vapor 
density  can  be  taken. 

For  the  range  of  20-80°C  the  data  can  be  approximated 
as  a  fourth  order  polynomial. 

Vv    -  (24  6   -   10.93    *    Tb    +    0.213    *   tb2   - 

0.00193    *    Tb3    +    6.55E-6    *    Tb4)    /   338 
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Tb     =       50°C    it    yields 
pv  =  10.67     (kg/m3 


where    Vv  -  — 

r  v 


5.  Specific  Heat  C0  (J/kg°C) 

For  all  Flourinert  Electrochemicals  (Figure  4) 
Cp    =    (0.241111    +    3.70374E-4    *    Tb)     *    4186 
at  T  =    50°C.  Calculation  yields  1086.7  J/kg°C. 

6.  Thermal  Diffusivity  a  (m2/s) 

a  =  k/pcp 
at  T  =  50°C,  calculation  yields  3.115E-8  m2/s. 

7.  Thermal  Expansion  Coefficient  (5  (1/°C) 
Table  4B  in  the  Product  Manual  shows 

P  =  0.00261    /    (1.740    -    0.00261)     *    Tb) 
at  T  =    50CC.  Calculation  results  in  0.00162  (1/°C. 

8.  Prandtl  Number  (Pr) 

Pr   =   v/cc 
at  T=50°C.This  yields  Pr  =  15.4 

9.  Calculation  of  Latent  Heat  hfg  (J/kg) 

Table  13A  of  the  Product  Manual  shows  the  values  can 
be  approximated  by  a  linear  equation. 

hfg   =    (8057   -   22.73    *    Tb)     *    10.79 
where  10.79  is  the  conversion  factor  from  Cal/Mole  to  (J/kg) . 
At  T  -  50°C  calculation  yields  74672  (J/kg) . 
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APPENDIX  B.  UNCERTAINTY  ANALYSIS 


UNCERTAINTY  IN  AREA 
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±  0.010  inches 
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-       Wire  diameter 

=   Length  of  heater  wire 
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3.04  *  10"5  m2 
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0.0001  \2  ,  /  0.020 
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CO, 


0.021    or    2.1% 


CO 


6.38    x    10"y    in'    =    9.89    x    10'4   m 


-4     ™2 


B.       UNCERTAINTY    IN    POWER 
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UNCERTAINTY    IN    HEAT    FLUX 
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D.   UNCERTAINTY  IN  TEMPERATURE 

Using  thermocouples: 
Variable 
Voltmeter  resolution 

Ice  reference  bath 

Temperature 

Polynomial  Temperature 

Conversion 

R 


Uncertainty 

0.026°C 

1.0    mv    resolution 

0.05°C 

0.00663°C   RMS 

1.0% 


cor. 


[  (C0W)2  +  (COIB)2  ]1/2  +CDC 


COT, 


[     (0.0025)2    +     (0.05)2    ]'1/2    +    0.00663 
0.063    or    6.3% 


E.   UNCERTAINTY  OF  WIRE  SURFACE  TEMPERATURE 

From  the  calibration  formula 


where, 
AT 


From  (1) 


R  =  R0   [  1  +  £  AT] 

Slope  of  the  calibration  curve 
Bulk  fluid  temperature. 


(1) 


R 
~R~ 


-  1 


/  e  +Tu 


For:  TR  =  44.14°C 
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TV    =    23.9°C 
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F.       UNCERTAINTY    IN   AT 

Where   AT  =    Ts    -    Tb 

dAT 
~6~T-  "        l 
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G.       UNCERTAINTY    IN    HEAT    TRANSFER   COEFFICIENT 
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q  is  the  heat  flux 
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APPENDIX  C.  CALIBRATION  OF  CHROMEL  WIRES 

The  chromel  wires  of  about  3  inch  length  were  immersed  in 
a  Rosemount  constant  temperature  distilled  water  bath,  and 
were  compared  with  a  platinum  resistance  thermometer  as  a 
standard.  The  temperature  was  raised  from  30°C  to  90°C  in 
increments  of  about  five  degrees.  At  each  increment,  the 
temperature  was  allowed  to  stabilize  for  five  minutes,  then 
measured  using  a  Newport  digital  thermistor  devise  for  the 
0.005  inch  diameter  chromel  wire.  All  recorded  data  are 
presented  in  the  following  table.  The  corresponding  best  fit 
line  is  also  presented  on  the  graph.  Two  sets  of  data  were 
taken  at  different  times. 


Thermistor 

T/C 

Resistance 

(C) 

(Ohm) 

29 

30.2 

4.74 

38.3 

39.7 

4.76 

49.5 

51 

4.78 

61.9 

63.4 

4.8 

76.2 

77.2 

4.83 

87.7 

88.4 

4.85 

52.3 

53.8 

4.79 

57.5 

59.1 

4.8 

60.5 
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4.8 
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4.81 
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72 
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4.82 

86.2 

86.9 

4.84 
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Figure  CI .  Chromel  Wire  Calibration 
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APPENDIX  D.  SOFTWARE 

Software  written  in  the  BASIC  language  converts  voltage 
measurements  into  temperature  values  for  the  bulk  fluid,  ice 
reference,  and  coolers  inside  the  thermocouples.  It  uses  the 
a  value  of  0.0004  Ohm/°C  to  calculate  the  surface  temperature 
offered  by  Kaye  (1976)  .  The  printed  output  covers  bulk  fluid, 
wire  surface  temperatures  and  heat  flux. 

A  FORTRAN  program  reads  the  output  values  created  by  the 
previous  program  and  calculates  liquid  properties  Ra,    Nu,    and 

Csf 

The  GRAPHER  program  graphs  the  data  and  prints  it  using  a 
HP-7440A  Plotter. 

Initial  experiment  data  was  fed  into  EASYPLOT  software  on 
the  Naval  Postgraduate  School's  main  frame  computer  via  floppy 
diskettes . 
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11 

OPTION  BASE  1 

12 

BEEF 

13 

PRINT  "INPUT  DATE" 

14 

INPUT  TanhS 

15 

PRINT  "  INFUT  TIME  " 

16 

INPUT  lam 

18 

(PRINT  "SCANNING  MEAS" 

19 

REAL  Erof(7  )  ,T<7>  ,0(3)  ,Uolt(  1  1  :19>  ,Pouier<  4  )  ,Hf  lux<  4  )  ,Ts(4>  ,0hm(4) 

2a 

REAL  Sumo 

25 

OUTPUT  70SrCGNFHEAS  DQU  ,10!  -107  ,USE  0' 

26 

FOR  1=1  TO  7 

27 

ENTER  709ttmf(I  ) 

- 

29 

i PRINT  ItEmf(I) 

30 

NEXT  I 

31 

OUTPUT  705 i "CONFHEAS  DCV  ,1 1 1-1 19, USE  0" 

32 

FOR  1=1 1  TO  19 

33 

ENTER  705; Volts  I  ) 

35 

i PRINT  I ;Uoit ( I  ) 

35 

NEXT  I 

37 

D<  i  )». 1 009605 i 

33 

D<2>=25727.9 

25 

D(3  )=-767345. 8 

40 

D(  4  ;=73025555 

41 

DCS  )=-9247455595 

42 

DC  5  )=6  .  95E-I- 1  i 

43 

0(7  )=-2.5SE+15 

44 

D(  5  )=3. 34E+ 14 

45 

FOR  1=1  TO  7 

46 

5uno=0. 

47 

FOR  J=1  TO  5 

45 

5uno=Sumo+D( j)«Emf(I)*(J-1  ) 

49 

NEXT  J 

50 

T(  I  )=5uno 

51 

Tb  =  (  T(  2 )+i  (  o  >+i ( 1  )  )/3 

52 

NEXT  I 

53 

PRINT  "DATE  ;";7arihS 

54 

PRINT  "TIME  ;  " ;  Zam 

55 

PRINT  "T" i i  ,T( i  ) ; °  C   UPPER  PORTION  OF  THE  FLUID" 

56 

(PRINT  USING  "T" ;1  ,DD.D(T( 1  );"  C   UPPER  PORTION  OF  THE  FLUID" 

57 

FRINT  "T";2  ,T(2);"  C   LOWER  PORTION  OF  THE  FLUID" 

53 

PRINT  •T,';3,T(3))"  C   LOWER  PORTION  OF  THE  FLUID" 

53 

PRINT  °T";4,T(4)i"  C   TE  COOLER  INSIDE  AL.  PLATE" 

66 

PRINT  "T";5,T(5>;"  C   TE  COOLER  INSIDE  AL.  PLATE" 

51 

PRINT  "T"|B,T(B)|"  C   AMBIENT  (ROOM  TEHPERATURE >" 

TEMPERATURE)" 
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55 
54 
55 
73 

as 

Si 
52 
S3 
54 

55 
5? 
98 
95 

isa 

1  Oi 
102 
i03 
i  04 
i  05 
1  57 
1  05 
1  03 
i  I  I 
1  12 
i  14 
1  15 
116 
I  17 
115 
1  19 
12© 
iZI 
123 
124 
1  25 
125 
!  27 
1  25 
12S 
1 33 
131 
132 
133 
1  34 
135 
i35 
340 


PRINT  "T" ; 7,7(7  ); "C    ICE  (REFERENCE   TEMPERATURE)" 

(NEXT  I 

FRINT  "TEULK-'iTb 

•power  coiculat ion3»»«»*»*»«*»» 

Fower<  1  )-< Volt(  12  )/.  1  >»Uolt(TD 

Power<3  >=■<  Vol  t  (  1 7  )/  .  1  >«tfoi  t  \  1"B  ') 

Iheater  surface  area  calculations 

!***•**♦••****•**•♦*#♦*♦•**»♦••••• 

Dia  = .005* .0254 

ni»3* .0254 

fiian-{ 3. 14-15325  >*Dia* A 1 

•heat  flux  calc. 

Hf  iux( 1  )»Power( 1  }/Aian 
Hf  lux(2  )=Foweri 2 >/Aian 
Hf iux( 5  5=Fower( 3  >/filan 
FRINT  "powerl ="  ;Fower( i  ; 
I  PRINT  "power2  =  '' ;  Fowerv  2  i 
SPRINT  "Dower3=" ;Power( 3  ) 

Ohm( I  >  =  Uoit( ll  >/< Voit<  «2  >/  .  1  ) 

!uhn(2 >=Voit( i  4  >/( voltv IS  )/. 1  ) 

Ohm  3  )=Uolii  IE )/( Voit (  17  )/' .  1  } 
; 0hn( 2  ) 
: 0hm( i  > 
;0hm(5  ) 


•FRINT  "0Hh'2  = 

•FRINT  "OHMi- 

! PR I NT  "OHM 3  = 

Dt=Tb-Ta 

1  Ta  15  THE  TEfir  OF  g'JLK  FLUID 


Tb  WHEN  POWER  OFF 


Ro  IS  THE  RESIS.  OF  CHRGnEL  WIRE  ACTUALLY  I 


7  TC 


I    THE   Ohm( 1  )    VALUE   WHEN    THE   SULK    FLUID    TEMP    ( Tb=Tamb  )    AND   FOWER   OFF 
I    (    POWER    OFF   MEANS    PowerC 1  >"0.0084    and    it    Is    obtained   by    setting    ) 
F0WER( 4 ) "0.0036 
R01=4.0ES  R04=5.793i 

(VIZ   40   V    Supplyist    to   00.1  0.i0   position   WHEN    1 5.4.    INSERTED 

Kaye    is    the    value    of    the   calibration    slope    RefiS.M    YOU   et    al] 

Ta  =  To 

Roi=4.QE= 

Ro3m3.7S81 

Kaye= . 0004 

R=4.055»(  l+,8004-*Dt  ; 

Ts(  1  >=( \ ( Onru i  )/Roi  )-l  )/Kaye  >+7a 

Ts<  3  )  =  (  (  (  Ohn<  3  )/F.o3  >- 1  i/Kaye  )  +  Ta 

Tf iln=( Tb+Ts( 4  )  )/2 

PRINT    "TSi=";Ts(  1  >,"  HFLUX    ;'-;Hflux(i) 

FRINT    "T53=";Ts(3)  ,"  HFLUX3    ; " ;Hf lux< 3  ) 

•liRINT    "T    filn    ;";Tfiin 

END 
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IHHHI  HHHMHHHHMMHHMHHMHHHH  IMMHMHMHHMHI  »HHHHHK«hhhhhh 

PPRPMETER  (  N=2S  ) 
INTEGER  I 


REAL  TS(N),TB(N),HF(N),NU(N),  RP(N),NUC(N)  ,  TF  (N)  ,  TF1  (N)  ,  BETP  (N) 
+      ,  DENF  (N)  ,  DENV  (N)  ,  KINV (N) ,  CP  (N)  ,  K  (N)  ,  HFG  (N)  ,  PLF  (N) ,  PR(N) 
+       ,SURT(N), RPYD(N), RPYFB <N) , NUSEX (N)  ,  NUSC (N) , DIP,  G,  V  (N) 
+       ,  MU  (N) ,  Y  (N) ,  X  (N) ,  US  (N) ,  TSPT,  CSF  (N) 

*  REPD  Tsurface.Heat  Flux,Tbulk  FROM  DPTP  FILE  • 

*  WRITE  Tsurface, Tbulk,Heat  Flux,  Nusselt, Ray leigh, Nusselt  (Chur)  » 

*  Tfilm  , beta, density, viscosity, Cp,  conduct ivity, latent  heat, 

*  Prant 1, surface  tension 

»  ♦^♦♦♦♦»»»»*HMMI^»H>H>'»«-^»»»-)H»***HHH»»*»»^  WHHHHHHHHMMHHW 

OPEN  (UNIT=1, FILE=*  INPUT. DPT'  ,  STPTUS='  OLD'  ) 

OPEN  (UNIT=2,  FILE='  OUTPUT.  DPT'  ,  STPTUS='  NEW  ) 

OPEN  (UNIT=3,  FILE='  PROPEL  DPT'  ,  STPTUS='  NEW  ) 

OPEN  <UNIT=4, FILE=' PR0PE2. DPT'  ,  STPTUS='  NEW  > 

OPEN  (UNIT=5, FILE=' ROHSE. DPT' , STPTUS=' NEW ) 

DO  15  I  =  1,N 

READ  (1,*)  TS(I),HF(I),TB(I) 
TF ( I )  =  ( TS  ( I > +TB ( I ) ) /2. 
TF1 (I)=TF(I)+273. 15 

BETR(I)=0. 30261/ ( 1. 740-0. 0026 1*TB ( I )  ) 
PRINT*, 'BETP  IS   =',BETP(I) 
DIP=0. 005*0. 0254 

G=9.ai 

TSPT=56. 

US  < I ) =TS ( I ) -TSPT 

DENF  ( I )  =  <  1.  740-0.  13026  1*TS  ( I )  )  *1000. 

V ( I )  = (246-10. 93*TB( I ) +0. 213*TB ( I ) **2-19. 3E-4*TB ( I ) **3+ 
+   6.55E-6*T8(I)**4) /338. 
DENV(I)=1. /V (I) 

KINV  ( I )  =  (0.  54-0.  01 16*TB  ( I )  +0.  i3002085*TB  ( I )  **2)  *lE-6 
CP  < I )  =  (0.  241 1 11+3. 7037E-4*TB( I) ) *4186. 
K  < I )  =  (-0. 001 15*TF ( I ) +0. 6033) /10. 
HFG (I)  =  (8057-22. 73*TB(I) )*10.  79 

PLF(I)=K(I)/(DENF(I)*CP(I) ) 
PR(I)=KINV(I) /PLF(I) 
SURT(I)=(-0. 12*TF(I)+15)/1000. 
MU ( I ) =K INV ( I ) *DENV ( I ) 

C  **********»HMt***-*************************************** 

Y(I)  =  (HF(I)/  (MU(I)*HFG(I)  )  )  *SQRT  (SURT  ( I )  /  (G*  (DENF  ( I) -DENV  ( I)  )  )  ) 
X  ( I )  =CP  ( I) * (TS ( I > -TB ( I ) ) / (HFG ( I ) *PR ( I ) **1 . 7) 
CSF ( I) =X ( I ) / (Y ( I ) **0. 333) 

C  ***-)M*****»HM»************HM»****************Ht»*»*************** 

*  CPLCULPTE  RPYLEIGH  and  NUSSELT  NUMBERS 

C        ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦■it-******************  »  »♦♦♦♦♦♦♦♦♦*«  *  »  »  ****** 
RPYD ( I ) =G*BETP ( I ) * ( TS ( I ) -TB ( I ) ) *D IP**3/ ( K INV ( I ) *PLF ( I )  > 
RPYFB ( I ) =G*BETP  ( I ) *HF  < I ) *DIP**4/ (K ( I ) *KINV ( I ) *PLF ( I) > 
NUSEX(I)=HF(I)*DIP/(K(I)*(TS(I)-TB(I)  ) ) 

NUSC ( I )  =  ( 0.  60+ ( 0. 3S7*RPYD ( I ) ** ( 1 .  /6.  ) / ( ( 1  +  ( 0. 559 /PR ( I )  >  **  ( 
+     9. /16.  )  )**(8. /27.  )  )  )  )**2 

WRITE  (2,20)  TS(I),TB(I),HF(I),NUSEX(I),  RPYD  ( I ) ,  NUSC  ( I ) 
WRITE  (3,30)  TF(I),BETP(I),DENF(I),DENV(I),KINV(I),CP(I) 
WRITE  (4,40)  K(I),HFG(I),PLF(I),PR(I),SURT(I) 
WRITE  (5,50)  Y(I), X(I),CSF(I) 
20       FORMRT  ( 1 X,  F8.  2,  2X,  F9.  1 ,  2X,  F12.  2,  2X,  F8.  3,  2X,  F12.  7,  2X,  F15.  3  ) 
30       FORMPT  (1X,F8.4,2X,F9.  4,2X,F9.  1,2X,F8.4,  2X,E10-3,2X,F15.  3  ) 
40       FORMPT  (1X,F8.  4,  2X,F14.4,  2X,E10.  3,  2X,F14.  4,  2X,F9.5  ) 
50       FORMPT  (1X,F8.4,  2X,F14.4,2X,E10.  3  ) 
15       CONTINUE 
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Theeis 

E56C5   Eren 

c.l       Heat  transfer  enhance- 
ment due  to  bubble  pump- 
ing in  FC-72  near  the 
saturation  temperature. 


